Magnetoresistance measurements of Graphene at the Charge Neutrality Point 
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We report on transport measurements of the insulating state that forms at the charge neutrality 
point of graphene in a magnetic field. Using both conventional two-terminal measurements, sensitive 
to bulk and edge conductance, and Corbino measurements, sensitive only to the bulk conductance, 
we observed a vanishing conductance with increasing magnetic fields. By examining the resistance 
changes of this insulating state with varying perpendicular and in-plane fields, we probe the spin- 
active components of the excitations in total fields of up to 45 Tesla. Our results indicate that v — 
quantum Hall state in single layer graphene is not spin polarized. 
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Under a magnetic field, the linear dispersion relation 
of low energy electron spectrum in graphene leads to 
unique Landau levels (LLs) whose energy difference is 
unequally spaced [IH3]. The LL spectrum, given by 
E n = ±^2nhvF 2 eB/c, where vp is the Fermi velocity 
and n = 0, ±1, ±2, ... is LL index, contains an n = 
level, termed the zero-energy LL (ZLL). In the absence 
of appreciable interactions or Zeeman splitting, each LL 
has a 4-fold degeneracy arising from a real spin and val- 
ley degeneracy. The appearance of the quantum Hall 
(QH) effect in graphene at the LL filling fractions v = 
±2, ±6, ... is a manifestation of this 4-fold degeneracy of 
graphene LLs [HE]. In the high magnetic field regime, 
however, this effective SU(4) spin-pseudospin symmetry 
can be broken, with more QH plateaus appearing at 
v = 0, ±1,±4 and developing signatures of QH states 
for other integer filling fractions [7-9 . The v = fill- 
ing factor that appears at the center of the ZLL presents 
something of a paradox in QH physics, as it is not marked 
by the usual longitudinal resistance minima that typify 
all other filling factors. While initial measurements on 
disordered samples at this filling factor reported high- 
field (above 30 T) resistance in the regime of tens of 
K^s [19 , subsequent reports on this quantum Hall state 
have shown a strong insulating behavior as sample mobil- 
ity is increased [TQHT3] , with two-terminal measurements 
of the highest mobility suspended samples measuring into 
the GfL range at fields as low as 5 Tesla [T4] H5], 

Theoretically, various models of symmetry breaking 
and ordering underlying this v = insulating state 
have been proposed. Most of the models fall under 
the framework of exchange-driven quantum Hall ferro- 
magnetism that separates different sectors of the SU(4) 
spin-pseudospin space [I6j [17]. These include: a fully 
spin-polarized ferromagnet [20] [23], a fully pseudospin- 
polarized charge density wave [22] [24], a Kekule dis- 
tortion with a spontaneous ordering of pseudospin [25] 
[26l [29] , and a canted antiferromagnet [27 : . An alterna- 
tive approach is based on magnetic catalysis: long-range 
electron-electron interactions that induce an excitonic 
gap [28] . Experimental reports on the nonzero filling frac- 
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FIG. 1: (a) Longitudinal conductivity (a xx in black) and hall 
conductivity(cr xy in red) as a function of back gate voltage 
at T = 1.7 K with B = 4 T; dashed arrows indicate QH 
states v — 0, 1, and 2; (b) Conductance (S) on a logarith- 
mic scale versus V g at B± = 3 T for two different tilting 
angles: 6 =18.3° and 6 =52.8°; (c) Maximum two-terminal 
resistance as B± increases for different tilting angles: 6 =18.3° 
(open triangle) , 61.8° (solid triangle), 72.2° (open square), and 
80.8° (solid square); the experimental error of is 0.3°; lower 
right inset: AFM image of the measured suspended graphene 
device; upper inset shows a schematic diagram for tilting an- 
gle and perpendicular and total magnetic field components. 



tion [8] suggest that the excitations of the v — 1 state 
have no spin, while the Kosterlitz-Thouless insulating be- 
havior of v = [11] is consistent with a Kekule distortion 
origin. The various models of the broken symmetry states 
involve unique bulk spin/pseudospin textures and corre- 
sponding edge state configurations [20] [21] . Thus trans- 
port measurements require careful comparison of the bulk 
and edge state conduction in order to answer questions 
related to the nature of the symmetry breaking at v = 0. 

In this letter we investigate the spin response of the 
v = QH state in monolayer graphene by measuring the 
bulk and edge conduction as a function of in-plane mag- 
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netic field using high- mobility suspended graphene and 
on-substrate graphene Corbino device. Our experiments 
reveal a vanishing conductance at v = 0, but neither ex- 
hibits an increasing gap with increasing in-plane field, 
suggesting that the v = state is not spin-polarized. 

The suspended graphene devices are prepared using 
the methods described in reference [30] : after thermally 
evaporating Cr/Au electrical contacts onto the mechan- 
ically exfoliated graphene samples [5], a chemical etch of 
buffered hydrofluoric acid is performed to remove the 
Si02 under the graphene sample, leaving the whole de- 
vice suspended approximately 200 nm above the SiC^/Si 
substrate. An atomic force microscope (AFM) image of 
the device is shown in the inset of Fig. 1(c). DC current 
annealing is then performed at low temperature (T = 
1.7 K) to remove residual impurities from the suspended 
graphene. Four-terminal transport measurements are 
conducted using conventional low-frequency lock-in tech- 
niques. The carrier density of the graphene is tuned by 
applying back gate voltage V g to the degenerately-doped 
Si substrate, with the magnitude of the tuned density 
determined using Hall measurements. The mobility of 
this annealed device is ~80,000 cm 2 /V-s. In Fig. 1(a), 
we show the longitudinal conductivity <j xx and Hall con- 
ductivity <j xy versus back gate at B =4 T normal to the 
graphene basal plane. As indicated by the vertical ar- 
rows, along with clearly developed v = 2 QH state, strong 
v = and developing v = 1 are observed as plateaus in 
a xy and the suppression of a xx at the corresponding fill- 
ing fractions. The appearance of the v — and v — ±1 
QH states indicates that the four-fold degeneracy of the 
ZLL is completely broken. 

To discern whether the v = symmetry breaking 
is spin-active, we apply a sequence of tilted magnetic 
fields that fix the perpendicular magnetic field B± while 
varying the total magnetic field B t . By fixing B± the 
magnetic length Ib = \/h/eB±, Coulomb energy scale 
E e - e — e 2 /Ait€Q€ r lB are held constant, meaning the 
electron-electron and exchange interactions that underlie 
the v = state are unchanged. However, if this state 
is fully spin polarized, the current-carrying excitations 
will have net spins that will be affected by changes in 
B t via the Zeeman energy AE Z = g{iBB t , where g is g- 
factor of electron and \±b is the Bohr magneton. At a 
fixed temperature the changes to the carrier excitation 
energy will result in a change in conductance observed 
at the v — filling factor. Thus by tuning only the Zee- 
man energy and examining changes in the conductance, 
we can determine if the activation of the v = state is 
spin-sensitive. 

The results of measuring the insulating state of the 
suspended device at several different tilting angles are 
shown in Fig. 1(c), where the resistance maximum Rmax-, 
is measured at the charge neutrality point V g = Vd, at 
a fixed base temperature T =1.6 K. Since the resistance 
for v = QH state tends to increase rapidly as a function 



of B in [TQl [TTJ [14j [15] , Rmax is a good measure to probe 
this insulating state. Here we use two-terminal current 
measurement with a constant voltage bias in order to 
eliminate any self- heating effects (< pW) and to maxi- 
mize the measurable resistance range. At T = 1.7 K, we 
found that R m ax increases from ~10 ¥St up to 100 MSI 
(comparable to the limit of our measurement set-up) as 
Bj_ changes from to 3 T. The tilting angle dependence 
of Rmax versus B± curves show such a trend: while we do 
not observe appreciable dependence of R m ax on in-plane 
magnetic field at lower values of the tilting angle (i.e., 
larger B±/B t ratio), there is an indication that Rmax 
decreases at larger 6 (i.e, smaller B±/B t ratio). This 
trend becomes most obvious for the largest tilting angle 
we measured, =80.8°, corresponding to B±/B t , where 
we observe that Rmax versus B± curve is substantially 
lower than any other curves in the graph. The observed 
trend in the suspended device, i.e., decreasing Rmax with 
decreasing B±/B t at fixed Bj_ suggests that the v = 
gap decreases as B t increases. This dependence can be 
viewed as strong evidence against a fully spin-polarized 
ordering of the v = QH state, as this ordering would 
result in an increase in the gap as B t increases. The rel- 
ative insensitivity of Rmax to changes in angle for the 
small tilt angles may be due to broadening induced by 
thermal smearing or disorders. 

There are two obstacles in using suspended samples to 
draw more quantitative conclusions about the nature of 
the v = QH state. First, due to the mechanical in- 
stability of suspended samples, Rmax drifts slightly with 
respect to V g . Fig. 1(b) shows the conductance as a func- 
tion of V g measured at two different tilting angles. Al- 
though the overall behavior is consistent, the position of 
V g where Rmax occurs is slightly shifted. Even worse, this 
shift changes when the device is thermally cycled, making 
it difficult to estimate the energy gap by the thermally- 
activated behavior. Second, the four-/two-terminal de- 
vice geometry measures both the bulk conductance and 
any possible edge conductance in parallel. This becomes 
a major source of ambiguity in distinguishing whether the 
observed insulating behavior originates from the bulk in- 
sulating state without the edge conduction or from the 
localization of edge states by spin/pseudospin-flip scat- 
tering [19, 20 . In order to avoid the mechanical instabil- 
ity and to isolate the bulk conductance, we employ an on- 
substrate Corbino geometry, a disk-shaped sample with 
coaxial contacts in which the current flows radially from 
an inner contact to an outer ring contact. This geometry 
not only eliminates any unknown edge effects that might 
interfere with determining the v = conductance, but is 
also insensitive to the formation of the known quantized 
edge conductances of other filling factors. This geometry 
then directly allows probing bulk conduction, and thus 
puts the v = insulator on an even footing with the bulk 
insulating character of every other filling factor [31]. 

The fabrication procedure for our Corbino devices is 
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FIG. 2: (a) Schematic side view of a Corbino device; (b) 
Optical device image before the inner electrode is contacted 
by the ground plate; (c) Optical image of a finished device; 
(d) Bulk Conductivity a xx as a function of back gate voltage 
V g at zero field (in red) and 14 T (in blue), at T =7 K. 
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FIG. 3: (a) <J xx cts a function of filling factor at four per- 
pendicular magnetic fields at T — 1.7 K. The shaded bands 
highlight the developing filling factors as the four-fold degen- 
eracy of the ZLL is broken, each manifested as a vanishing 
bulk conductance. 



shown as in Fig. 2(a). Monolayer graphene pieces are de- 
posited on SiO2(300nm)/Si substrates using established 
mechanical exfoliation techniques, then Au/Cr ring-like 
electrodes are fabricated by e-beam lithography, (an op- 
tical image is shown in Fig. 2(b)), followed by a dielectric 
layer deposition and a top Au/Cr plate contact to con- 
nect to the inner contact (as shown in Fig. 2(c)). The 
plate geometry connecting to the inner contact guaran- 
tees that any voltage applied to this contact will result 
in a uniform change to the graphene carrier density. To 
measure the bulk conductance of the graphene, we ap- 
ply an AC voltage bias (Vuas) across the inner and outer 
contacts, and measure the current (J) using a current 
preamplifier and lock-in amplifier. The bulk conductivity 
is then given by a xx = (\n(r out /r in )/27r)(I /V bias ), where 
r ou t and Ti n are the radii of the outer and inner contacts, 
respectively. 

Changing the back gate voltage V g: we can tune the 
carrier density in the graphene channel connecting the 
inner and outer contacts of the Corbino device. Fig. 2(d) 
shows the bulk conductivity <j xx vs. back gate voltage 
V gi at B =0 T and 14 T at temperatures lower than 
7 K. The mobility of this particular Corbino device is 
^13,000 cm 2 /V-s, obtained from the zero- field resistance. 
At B = 14 T, the four-fold degenerate QH state filling 
factors v = ±2, ±6, ±10 appear as vanishing cr xx at their 
corresponding carrier density. The gate capacitance of 
this device is estimated to be C g /e = 7.1 x 10 10 cm _2 V _1 
from the position of the observed conductivity minima. 

Since the mobility of the on-substrate Corbino devices 
is lower than that of the suspended devices, relatively 
higher magnetic fields are required to access the degen- 
erately broken filling factors. As shown in Fig. 3, at 



low field (B = 11.5 T) well-defined v — ±2 states are 
observed on both sides of the charge neutrality point, in- 
dicative of the four-fold QH degeneracy. As the magnetic 
field increases to 18 T, a dip of bulk conductivity appears 
at the charge neutrality point. This dip fully evolves and 
the current flow falls below the noise level at B = 30 T. 
This observation of a vanishing bulk conductivity is con- 
sistent with the formation of the v = QH state [31] . At 
the same magnetic field, the conductivity minima corre- 
sponding to the v = ±1 filling factors are visible. At 
B = 45 T, the four- fold degeneracy at the zero energy 
level is completely lifted, and the LL splitting at v = —4 
that marks the degeneracy breaking of the n = 1 LL is 
apparent, similar to the previous observation [7]. In all 
measured devices, the magneto-conductance is strongly 
suppressed in the regime between the v = — 1 and v — — 2 
filling factors and is not measurable within our experi- 
mental sensitivity, which remains not fully understood. 

As with the suspended devices, we adjust the rela- 
tive strengths of the Zeeman and Coulomb energy in the 
Corbino devices by tilting the field in order to explore the 
nature of the v = degeneracy breaking. In Fig. 4(a), 
<j xx vs. filling factor v is plotted with constant normal 
field (Bj_ = 21 T) and with the total field (B t =Bj_/ cos 0) 
increasing. Taking the dielectric constant e r =4, the char- 
acteristic Coulomb interaction energy at B± = 21 T is 
E e - e = 740 K, while the Zeeman energy varies from 
E z = 47 K at B t = 35 T to E z = 60 K at B t = 45 T. As 
the Zeeman energy is increased, the behavior of the v = 
and v = 4 states are completely different. For the v — ±4 
QH state, the <j xx minima decrease with increasing B tl 
indicating that a spin polarization underlies this LL, a 
finding consistent with previous experiments on Hall bar 
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FIG. 4: (a) a xx as a function of filling factor at a constant 
perpendicular magnetic field B± = 21 T for several different 
total magnetic field. Data taken at T — 1.7 K. From top to 
bottom the corresponding total fields are Bt = 35, 40, and 
45 T, respectively. The inset shows a close-up of the main 
panel at the v — QH state, (b) <j xx at B± = 14 and 15 T, 
where it is sensitive to small changes in the perpendicular 
field. 

devices [7]. In contrast, the conductance curves of the 
v = state coincide with each other as the total field is 
increased from B t = 35 T to 45 T. The conductance min- 
ima are unvarying even in a magnified logarithmic-scale 
view, as shown in the middle inset of Fig. 4(a). The fact 
that this minima is independent, within disorder broad- 
ening, to changes in the in-plane field is also consistent 
with a state that is not fully spin-polarized, and adds fur- 
ther credence to the hypothesis that the v — symmetry 
breaking is not of spin origin. 

We also perform fine-tuned tilted field measurements 
in a range where the change of Zeeman energy is larger 
(increased by 50%) and the v = minima are more sen- 
sitive to small changes in B±. Fig. 4(b) shows a log-scale 
<j xx vs. filling factor v at B± = 14 T and 15 T. As the 
normal field increases by ~ 6%, there is a decrease in 
the bulk conductivity minima, showing that the v = 
state is not yet fully developed. Increasing the total field 
by ~ 50% while fixing B±, the minima display the same 
insensitivity to in-plane field as in Fig. 4(a), reaffirming 
that the excitations of the v = state have no net spin. 
We are aware of not observing a decreasing R m ax with de- 
creasing B±/B t at fixed B± in Corbino device, although 
the range of the Zeeman to Coulomb energy ratio change 
in Corbino device is similar to that of the suspended de- 
vice. The discrepancy of the behaviors could be under- 
stood to be consequences of the different disorder energy 
scale. 

As to the v = 1 QH state, experimental data of its 
tilted- field dependence is also shown in Fig. 4(a). The 
<j xx minima at this filling factor decrease as B t increases. 



This observation implies that the origin of this state is 
due in part to a lifting of real spin degeneracy. Combining 
this observation with that of the v = state, it produces 
a symmetry-breaking picture of the ZLL where a non- 
spin polarized state forms at v = and a spin-polarized 
state with spin- flip excitations forms at v — 1 [32] . 

We are also aware that the observation of the spin- 
active v = 1 character is inconsistent with the obser- 
vations of Jiang et al. [8 , whose measurements implied 
that excitations of v = ±1 QH states has not spin flip. 
This raises the possibility that the excitations at v = 1 
and its ground state may depend on the specific disorder 
concentration in individual samples p~6j [18]. However, 
the fact that the insulating v = state does not respond 
to increasing in-plane magnetic fields in both suspended 
and Corbino devices, where disorder densities are very 
different, provides evidence that disorder effects do not 
alter our conclusion that the v — QH state is not a 
spin-polarized state for a wide range of disorder. 
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